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ABSTRACT 

Microwave  radiation  is  absorbed  by  atmospheric  oxygen  and  water 
vapor  and  is  absorbed  and  scattered  by  rainfall.  This  report  provides 
values  of  these  losses  in  the  frequency  decade  from  10  to  100  GHz  for 
paths  through  the  atmosphere  at  elevation  angles  from  0  to  90  degrees. 

The  theoretical  express*ons  for  oxygen  absorption  arc  assembled 
in  explicit  form.  The  resonant  frequencies  are  taken  from  published 
low-pressure  measurements,  and  the  entire  absorption  profile  is  fitted 
to  published  intermediate  and  atmospheric  pressure  results,  with  a 
constant  of  proportionality  and  the  absorption  resonance  line  width  used 
as  fitting  parameters.  The  absorption  line  breadth  is  found  to  deviate 
from  the  accepted  linearpressuredepender.ee,  increasing  approximately 
as  the  square  root  of  pressure  for  pressures  above  one-ha  if  atmospheric. 

The  results  of  similar  analyses  of  water  vapor  absorption  are  as¬ 
sembled  and  combined  with  the  results  of  oxygen  absorption  analyses 
to  obtain  the  complete  absorption  coefficients  as  a  function  of  the  at¬ 
mospheric  parameters. 

The  total  attenuation  aue  to  absorption  in  the  clear  atmosphere  is 
computed  for  a  model  atmosphere  with  various  water  vapor  concentra¬ 
tions.  Results  are  provided  in  graphical  form,  giving  attenuation  as  a 
function  of  frequency  and  elevation  angle. 

Tiie  pre-diction  of  rainfall  attenuation  is  severely  limited  by  the 
absence  of  meaningful  statistics  on  rain  rate  and  storm  structure. 
Rough  estimates  for  rain  attenuation  are  made  for  crude  model  storms 
which  occur  for  1,  0.1,  and  0.01  percent  of  the  time.  Total  values  of  the 
attenuation  for  the'  model  storms  and  the  atmospheric  absorntion  are 
provided. 


PROBLEM  STATUS 

This  is  an  interim  report  on  one  pliase  of  the  pi  oblems;  work  on 
other  phases  is  continuing. 
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ATTENUATION  OF  MICROWAVE  RADIATION  FOR 
PATHS  THROUGH  THE  ATMOSPHERE 


INTRODUCTION 

The  attenuation  of  microwave  radiation  in  the  frequency  dtcade  10  to  100  GHz  for 
propagation  paths  through  the  atmosphere  at  all  elevation  angles  is  of  great  interest  in 
the  analysis  of  experiments  in  radio  astronomy  and  communications  by  satellite.  The 
purpose  of  this  report  is  to  assemble  the  explicit  theoretical  expressions  and  measured 
constants  used  to  compute  the  attenuation  and  to  orovide  the  results  of  such  a  calculation 
for  a  typical  model  atmosphere.  The  computation  of  attenuation  for  paths  through  the 
atmosphere  requires  a  knowledge  of  the  magnitude  of  the  absorption  by  the  atmospheric 
constituents,  the  frequency  dependence,  and  the  spatial  variation  of  the  effect  along  the 
path.  The  attenuation,  in  general,  consists  of  two  contributions  that  require  quite  different 
treatments.  The  first  contribution  is  the  absorption  by  the  molecular  oxygen  and  water 
vapor  that  controls  the  attenuation  in  the  clear  atmosphere  and  can  be  computed  accu¬ 
rately.  The  second  contribution  is  attenuation  due  to  rainfall;  this  can  be  treated  only 
qualitatively,  but  it  is  quite  important  to  the  communicator  who  is  not  always  able  to 
choose  ideal  conditions. 

Molecular  oxygen  and  water  vapor  absorption  resonances  or  lines  completely  domi¬ 
nate  the  absorption  spectrum  in  the  band  from  10  to  100  GHz.  A  review  of  the  existing 
theoretical  treatments  and  experimental  measurements  is  made  here  to  obtain  expres¬ 
sions  giving  the  coefficient  of  absorption  as  a  function  of  frequency,  temperature,  and 
pressure.  The  total  attenuation  for  paths  of  varying  elevation  angles  is  obtained  for  a 
model  atmosphere  by  integrating  the  semiempirical  expressions  along  the  path. 

First,  the  theoretical  analyses  are  used  to  obtain  the  general  expression  describing 
the  interaction  of  the  incident  radiation  with  the  oxygen  molecule.  Then  the  various 
terms  in  the  expression  are  identified  and  reduced  to  the  form  to  be  used  in  the  path 
computations.  In  its  final  form,  the  expression  determines  the  functional  dependence  of 
the  absorption  spectrum  except  for  a  function  describing  the  shape  of  the  resonances  in 
the  frequency  domain  and  a  constant. 

The  expression  obtained  for  oxygen  is  adapted  to  the  essentially  similar  water  vapor 
interaction,  and  a  semiempirical  approximation  of  the  contribution  of  resonances  that 
occur  beyond  lOu  GHz  is  introduced.  The  problem  of  obtaining  a  function  to  describe  the 
shape  of  the  resonance  line  is  discussed,  introducing  additional  pressure  and  temperature 
dependence  and  another  fitting  parameter,  the  line  breadth. 

The  constants  which  remain  in  the  theoretical  functions  are  obtained  by  fitting  pub¬ 
lished  experimental  data.  Evidence  for  a  pressure  dependence  for  the  oxygen  line 
breadth  different  from  the  theoretical  result  is  also  derived  5m  the  data.  The  total 
attenuation  due  to  absorption  by  oxygen  and  water  vapor  is  "->uted  as  a  function  of  fre¬ 
quency  and  elevation  angle  for  the  model  atmosphere. 

The  attenuation  due  to  rainfall  is  considered  in  an  approximate  manner  in  the  last 
part  of  this  report.  The  theoretical  values  of  the  attenuation  as  a  function  cf  frequency 
and  rain  rate  are  taken  from  published  results.  The  severe  limitations  of  the  existing 
statistical  data  on  the  spatial  and  temporal  distributions  of  rain  rate  are  discussed,  and 
some  crude  model  storms  are  constructed.  Finally,  the  attenuations  produced  by  the 
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rain  for  propagation  paths  through  the  model  storms  are  computed  and  combined  with  tne 
attenuation  due  to  oxygen  and  water  vapor  to  obtain  representative  values  of  the  total 
losses. 


ABSORPTION  BY  MOLECULAR  OXYGEN 

The  absorption  of  microwave  energy  by  oxygen  results  from  magnetic  dipole  inter¬ 
actions  with  the  incident  radiation  due  to  the  oxygen  molecule's  permanent  magnetic 
dipole  moment,  which  produces  transitions  between  molecular  fine  structure  levels  ui 
the  allowed  rotational  states.  A  detailed  discussion  of  the  theory  of  microwave  spec¬ 
troscopy  is  beyond  the  scope  of  this  report,  and  the  various  treatises  (1-5)  on  which  this 
work  is  based  should  be  referred  to  for  an  exact  exposition. 

The  actual  amount  of  the  incident  radiation  absorbed  by  a  given  transition  in  a  gas  of 
the  molecules  depends  on  the  relative  fraction  (/, )  of  the  molecules  which  are  in  the 
lower  energy  state  of  the  transition  and  on  the  square  of  the  matrix  element  of  the  dipole 
moment  connecting  the  two  states.  The  magnitude  of  the  matrix  element  is  a  measure  of 
the  effectiveness  of  the  radiation  in  producing  the  transition. 

The  complete  expression  as  obtained  by  Van  Vieck  (4)  for  the  absorption  coefficient 
y  of  a  gas  of  molecules  having  a  series  of  transitions  has  the  form 


y  -  (constant)  f N  ^ 
'  .  ; 


f,S(v,x-y) 


I' 


This  expression  consists  of  pairs  of  terms  representing  emission  anti  absorption  transi¬ 
tions  between  two  states.  The  algebraic  sum  of  the  two  terms  is  the  net  absorption  and 
is  dependent  on  the  relative  population  of  states  of  higher  and  lower  energy.  The  sub¬ 
scripts  t  and  ;  are  labels  representing  states  defined  by  a  set  of  quantum  numbers.  The 
frequency  of  the  incident  radiation  is  v,  and  v  t  is  the  resonant  frequency  of  the  absorp¬ 
tion  line.  The  resonant  frequency  is  the  frequency  of  the  quanta  of  radiation  energy 
equal  to  the  difference  between  the  energies  of  the  initial  and  final  states  of  the  transition. 
N  is  the  number  of  molecules  per  unit  volume,  s  ( i/; , .  u)  is  a  function  describing  the  ab¬ 
sorption  line  profile  (i.e.,  the  absorption  versus  frequency  curve),  is  the  matrix  ele¬ 
ment  connecting  states  i  and  / .  The  constant  is  a  combination  of  physical  and  mathe¬ 
matical  constants  and  is  given  by  the  exact  theoretical  treatments.  As  the  predicted 
value  does  not  fit  the  data  perfectly,  the  constant  will  be  used  as  a  fitting  parameter,  and 
any  physical  constants  subsequently  introduced  in  this  discussion  will  be  absorbed  into 
this  term. 

The  energy  level  structure  of  the  oxygen  molecule  with  intervals  in  the  microwave 
region  is  produced  by  interaction  of  the  uncompensated  electronic  spin  (total  spin  for 
molecule  =  1)  and  the  "end  over  end"  rotational  angular  momentum  of  the  molecule.  The 
rotational  levels  (quantum  number  K )  are  split  into  levels  corresponding  to  the  allow¬ 
able  ways  of  constructing  the  total  angular  momentum  vector.  This  level  structure  is 
modified  further  by  the  rotation  of  the  molecule,  the  interaction  between  the  processing 
component  of  the  electron  orbital  angular  momentum,  and  the  total  spin  and  centrifugal 
deformation  of  the  molecule, 

Symmetiv  considerations  allow  the  rotational  quantum  number  A  to  take  on  odd  in¬ 
tegral  values  1.3,...,  only.  Corresponding  to  each  A',  two  possible  transitions, 

J  -  K  -  i  -  J  -  A'  and  J  -■  K  -  J  -  A'  *  l  (J  =  total  angular  momentum  quantum  number),  are 
allowed  by  the  selection  rules  A/  =  ♦  i ,  AA'  =  o  .  The  frequency  of  the  quanta  absorbed  or 
emitted  can  then  be  calculated  from  the  difference  in  the  energy  levels.  The  fine  struc¬ 
ture  of  oxygen  results  in  an  infinite  series  of  lines  grouped  about  60  GHz,  of  which 
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approximately  25  can  be  detected  in  actual  measurements.  The  energy  levels  have  been 
computed  by  Schlapp  (3)  in  fair  detail,  and  Burkhaiter  et  al.  ',6)  have  obtained  the  fre¬ 
quencies  of  25  lines,  which  could  be  fitted  to  the  Schlapp  theory  with  the  addition  of  em¬ 
pirical  terms.  The  frequencies  of  additional  lines  too  weak  to  measure  can  be  obtained 
from  such  semiempiricai  fi's. 

The  matrix  elements  are  derived  through  quantum  vector  addition  and  obtain  the 
form 


“ K »’  f  2  A'  +  I ) !  .V  -  l> 

2  (  A  -  1)'2A-  I) 

'**  (  2  A'  •  I  >  A' 

2  8,.0‘  l  A2  *  A  f  I )  (  2A  *  1 ) 

"A„  A!  A'  -  1  M  2A  -  ll  ’ 

where  -0  is  the  Bohr  magneton.  It  should  be  noted  that  oxygen  has  lines  at  , )  =  o  for 
each  rotational  level,  which  result  in  a  nonresonant  contribution  to  >  total.  The  subscripts 
A',  and  A',  denote  the  transitions  I  -  h-i—J-K  and  ./  -  x-J  -  A  ?  l  respectively,  while 
A0  denotes  the  nonresonant  term  corresponding  to  that  value  of  A'. 

The  fraction  of  molecules  in  the  /th  state  can  be  expressed  as 

f .  -  U.f.i  If  e»P  (  -  A',  k  T)  . 

Z  =  sL>  (*.[.)[  c\ p  t -  A'2  kT)  . 


a  Boltzmann  distribution  with  statistical  weight  factors  (w.f.)  reflecting  the  more  fre¬ 
quent  occurrence  of  certain  states  due  to  degeneracies.  Z  is  the  partition  function,  k  is 
Boltzmann's  constant,  T  the  absolute  temperature,  and  K,  the  energy  of  the  ;th  state, 
which  depends  in  general  on  all  of  the  quantum  numbers  defining  the  state. 


In  detailed  treatments  (see,  for  example,  Ref.  2),  f ,  is  expressed  as  /,.<>> J  ( J.  A',  i/>, 

in  which  each  term  has  the  form  given  above;  /,  is  the  fraction  of  states  in  a  vibrational 
state  defined  by  the  principle  quantum  number  n  ,  and  f  r  the  fraction  of  states  which  are 
in  a  rotational  state  defined  by  A',  and  i/.  For  the  oxygen  molecule,  the  lowest  vibra¬ 
tional  frequency  is  sufficiently  high  that  fv  is  nearly  1  for  the  vibrational  ground  state  at 
temperatures  experienced  in  the  atmosphere,  and  f .  can  be  considered  equal  to  fr  . 

If  the  identity -  /,  /  ,<  i  -  ft  f  )  is  used,  then  the  population  difference  is 


\ 


-a  tr 

(w.f  I  f  ’ 

-t  \T 
(w.f  ’ 


For  transitions  between  states  of  constant  A  . 

U  |-.)t  =  (w  f.l  -  2A‘  1  . 


and 
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-(r.  )  lr  ■ 

'  J 


or 


and  since  A t  , ;  <<  *r, 


-  •  •  /, 

/  -  /  -  /  [i  - "Al  <  >  * rl 
;  '  ;  L  j 

A  1/ 


The  dependence  of  >  on  T  obtained  by  Van  Vleck  (4)  further  assumes  tne  energy  K  to 
depend  cn  K  only,  and  thus  f  has  the  approximate  form 


/;  *  h 


< : -v *  nr 


-  /■'.  » r 


y  kT 

L.  f  2/*  !)<•  • 


where 

kk  =  v.k  ( A  .  n  and 
H  =  2.862  x  10- 21  joule. 

Since  lri;!l5  =  ;  ar.d  v>  =  -$<  ✓,,.:)  (see  Eq.  (3)),  >  car.  oe  written 

t  -  ( constant ) vK  L 

t  , ; 

and  y  becomes 

y-  (roootact)  /  !m,;!  v,  ,  St  ^ ,  ,  •  *• )  f, 

■ 


(1) 


If  the  atmosphere  is  assumed  to  be  sufficiently  dilute  to  be  approximated  by  an  ideal  gas 
v  =  /'0  kT  where  P0  is  partial  pressure  of  oxygen,  so  that 


y  =  (  const.ioi  ) 


« ; 


\x  the  pressure  dependence  introduced  through  S(*t}.v)  is  neglected  for  the  moment,  > 
is  directly  proportional  to  partial  pressure.  At  moderate  temperatures,  the  energy  lev¬ 
els  are  small  compared  to  kT.  E  %  o  03  kT,  and  the  sum  over  all  states  in  the  denomi¬ 
nator  of  ft  is  approximately  proportional  to  r.  At  frequencies  removed  from  the  reso¬ 
nant  frequencies,  the  temperature  dependence  is  given  by  7” 2  because  many  lines 
contribute  to  the  absorption  so  that  the  sum  of  the  exponentials  in  the  numerator  tends 
to  cancel  that  in  the  denominator  (except  for  additional  dependence  introduced  through 
St  »<  .  ).  At  frequencies  where  one  line  dominates,  the  dependence  is  given  by 


1 1 


exp  (  -  E j>  kT) 


The  form  finally  used  in  computing  the  absorption  coefficients  is 
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><;•.  T,  lb 


/']  7ol  v.*l 

>o  »T  r  *  L  ,  121  '  "  *H 

L  °J  L  j  L/tr.,1 


v 

A  odd 


A  (  2  A  -*  ‘i )  r  i.,  >i 

rrr  KMiav^ 


(A'»  li(2A-  It 


-  [:'A  5 1 .  K  .  v)] 


2(  A-2  t  A  f  I )  (  IK  *  I )  < 


Hi  A'  *  I ) 


J  I  i  m  [ i *N<  x.  i')  j  1  ^  exp  ( -  kT) 

I  j  -  3  j  I 


(2) 


>0  is  one  of  the  constants  adjusted  in  this  work  to  give  the  best  fit  to  the  experimen¬ 
tal  data  at  l\,  and  7'0 1  the  ;  K  are  the  frequency  intervals  or  resonances  of  the  spectrum, 
and 


Kh  k-  =  (2  073  I)  a<A’  )  <  "Kt  . 

The  summation  is  over  odd  values  of  K  only  and  is  cut  off  at  A  =  51,  as  the  contribution 
of  lines  beyond  this  was  found  to  be  insignificant.:  The  >k.  are  taken  from  experimental 
results  discussed  in  a  later  section  for  A  =  1  to  25  and  from  the  semiempirical  models 
for  A'  =  25  to  51:  * 


uK  ( Gib)  =  .39.550  -  0.252A'  -  910  A 
(Gib)  -  39  790  *  0  2T>2A'  -  9  10  A' 


ABSORPTION  BY  WATER  VAPOR 


The  absorption  of  microwave  energy  by  water  vapor  (4,7)  results  from  electric 
dipole  interactions,  and  the  theoretical  treatment  takes  the  same  general  form  as  that 
given  for  oxygen  in  the  previous  section. 


The  rotational  level  structure  of  water  vapor  is  much  more  complex,  consisting  of 
a  multitude  of  levels  spaced  at  intervals  in  the  microwave  regions.  Fortunately,  the 
selection  rules  that  result  from  the  symmetry  properties  of  the  molecular  levels  permit 
only  one  transition  in  the  frequency  range  of  10  to  ICO  GHz  (at  :  0  =  22.254  GHz),  while 
several  quite  intense  lines  occur  in  the  far  infrared  (;  >  100  GHz).  The  absorption  by 
water  molecules  then  takes  the  form  of  one  term  of  the  sum  of  Eq.  (1)  above  plus  an  ad¬ 
ditional  semiempirical  term  approximating  the  total  contribution  of  the  low-frequency 
tails  of  the  higher  frequency  lines.  The  fraction  of  molecules  in  the  ;th  state,  /,  ,  has 
the  form  y.~'  exp  (-  ,  k n ,  where  /.  is  again  the  partition  function  which,  for  the  case  of 
water  vapor,  is  a  double  sum  over  two  quantum  numbers  with  significant  contributions 
from  many  states.  By  numerical  evaluation,  it  can  be  shown  that  7.  *  r'  -  within  a  10- 
percent  error  for  variations  of  about  120°K  in  temperature  (4). 

The  number  density  of  water  vapor  molecules  \  (s  proportional  to  the  mass  den¬ 
sity  .  .  The  water  vapor  absorption  can  be  expressed  in  terms  of  the  vapor  density  . 
and 


,.(-a;p  kD 


.  iO  . 


where  is  one  of  the  constants  adjusted  to  fit  the  experimental  results. 


♦  Adapted  from  the  wave  number  result  given  in  Table  J  of  Ref.  4. 
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The  residual  absorption  of  higher  frequency  lines  is  given  by  a  sum  over  the  lines 
with  replaced  by  its  l*t:it  at  y,,  >>  y.  The  exact  calculation  of  this  term  is 

unprofitable,  as  the  result  is  found  to  be  only  one-fourth  of  the  observed  residual  ab¬ 
sorption,  probably  due  to  the  inadequacy  of  the  theory  of  S(  >, , . :  >  in  this  limit.  The 
semiempirical  result  of  Van  Vleck  (5)  adeouately  describes  the  experimental  results  and 
will  be  used  for  the  present  purpose: 


ABSORPTION  LINE  SHAPE  FUNCTION 

The  frequency  dependence  of  the  absorption  coefficients  obtained  as  discussed  pre¬ 
viously  is  largely  determined  by  the  line  shape  function  . .  ■ )  .  The  natural  widths  of 
the  absorption  resonances  are  very  small  (~  10 Hz),  resulting  only  from  the  growth  of 
the  amplitude  of  the  vibrations  upon  absorption  of  the  incident  quanta.  The  observed  line 
widths  are  of  the  orde;  of  109  Hz  at  atmospheric  pressure,  and  the  absorption  at  fre¬ 
quencies  removed  from  the  resonant  frequencies  is  completely  controlled  by  this  line 
broadening. 

In  the  case  of  atmospheric  gases,  the  lines  are  broadened  by  collisions  between  the 
molecules.  Van  Vleck  and  V/eisskopf  (8)  have  treated  this  classically  by  assuming  the 
collisions  to  be  sufficiently  strong  that  the  oscillator  is  stopped  at  the  moment  of  colli¬ 
sion  and  later  restarted,  with  its  phase  with  respect  to  the  incident  field  determinec  by 
the  condition  that  thermodynamic  equilibrium  exists  between  the  ensemble  of  oscillators 
and  the  field.  This  requires  that  the  time  of  the  collision  is  much  less  than  the  mean 
time  between  collisions.  The  kinetic  theory  is  applied  to  obtain  the  probability  of  occur¬ 
rence  for  collisions.  The  results  are  combined  and  averaged  over  all  the  molecules  to 
obtain  the  susceptibility  and  thereby  the  imaginary  part  of  the  index  of  refraction,  the 
absorption  coefficient. 

Tneii  result,  known  as  the  Van  Vleck-Weisskopf  shape  factor,  is 


(3) 


where  Av  is  the  line  breadth  factor  A^  -  !  2-'  (-  is  mean  time  between  collisions),  which 
can  obtain  values  of  the  order  of  , . 


The  limit  used  to  describe  the  nonresonant  absorption  of  oxygen  is 


llltl  U'5l  -T.  A  J  -  — — — — - . 

.2  A  ,  2 

The  limit  of  S(>  .*)  used  to  describe  the  absorption  of  the  submillimeter  water  vapor 
lines  is 


( i  ons  t.int  )  ivA;/  . 

V3  *  >  *' 

Others  have  treated  the  problem  using  the  methods  of  quantum  mechanics  and  have 
obtained  the  same  result  when  the  same  assumptions  were  made.  There  is,  however, 
some  question  as  to  the  validity  of  the  simple  assumption  of  strong  collisions,  and  many 
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of  the  possible  intermolecular  forces,  some  acting  at  relatively  long  range,  may  affect 
the  profiles  of  the  absorption  lines.  In  experimental  determinations  of  -'V,  it  is  in  gen¬ 
eral  difficult  to  fit  the  entire  profile  with  this  line  shape  function  for  a  single  value  of  Ai/; 
for  example,  if  '  :  is  chosen  to  fit  the  region  contiguous  to  the  peak,  the  theoretical  ex¬ 
pression  does  not  fit  tne  tail. 

Since  the  available  empirical  values  of  actual  atmospheric  absorption  can  be  fitted 
to  within  their  uncertainties  by  appropriate  adjustment  of  the  constants,  the  Van  Vleck- 
Weisskopf  line  shape  is  usable  in  the  present  effort  as  a  good  model  even  though  the  the¬ 
oretical  treatment  may  not  be  completely  satisfying. 

The  line  breadth  factor  depends  on  the  total  number  of  molecules  per  unit  volume  N , 
their  velocity  e  ,  and  the  cross  section  for  collision  A.  - .  The  collision  cross 
section  can  be  shown  to  be  proportional  to  '  "* ' '  ,  where  n  is  the  exponent  of  the 
intermolecular  force  law  (P  *")  and  can  assume  values  from  1  to  7.  From  kinetic  the¬ 
ory,  v  is  proportional  to  7'1  • ,  and  from  the  ideal  gas  law  \  ■>  /*  T,  so  that 

n  •»  I 
2  <  r*  -  I  ) 

Or,  if  <  ‘  >  0  is  the  line  breadth  measured  at  P0 ,  lu ,  then 


where  r,  the  exponent  of  the  temperature  dependence  at  constant  pressure,  is 

it  -  i 
2<n  -  1) 

The  simple  linear  pressure  dependence  of  -V.  is  tested  on  experimental  data  in  the 
next  section  and  is  found  to  hold  only  at  lower  pressures.  Apparently  at  pressures  of 
300  mm  Hg  and  higher,  the  individual  collisions  become  less  effective  ir>  broadening  the 
line  (possibly  due  to  either  an  increasing  number  of  multiple  molecular  collisions  or  the 
failure  of  a  molecule  to  regain  equilibrium  between  collisions),  with  the  result  that 
measurements  in  the  range  of  300  to  750  mm  Hg  appear  to  exhibit  an  approximately  P'  3 
dependence. 

The  pressure  and  temperature  dependence  of  the  line  breadth  factor  enters  into  the 
complct'  expression  for  the  absorption  coefficient  (Eq.  (2))  through  s<: . . ,i ' ,  contribut¬ 
ing  a  further  difference  in  the  pressure  and  temperature  variation  of  ■(-..TP)  between 
that  at  frequencies  adjacent  to  the  resonant  peaks  and  that  at  frequencies  removed  from 
the  peaks.  This  simple  behavior  of  the  line  breadth  is  further  complicated  by  the  de¬ 
pendence  of  the  perturbing  effect  of  the  collisions  on  A  .  Measurements* on  individual 
lines  at  low  pressures  show  that  the  variation  of  effective  width  from  line  to  line  is  less 
dian  10  percent  of  the  width,  so  that  using  an  average  value  weighted  over  line  strength 
results  in  a  small  error  (3). 


DETERMINATION  OF  CONSTANTS  BASED  ON 
PUBLISHED  EXPERIMENTAL  DATA 

Oxygen 

The  experiments  that  have  beer,  conducted  are  of  two  general  types:  measurements 
of  ) ,  the  i-A  or  V  at  reduced  pressure  and  measurement  of  the  actual  effect  of  the 
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absorption  at  atmospheric  pressure  on  a  signal  —  either  in  the  laboratory  or  over  a 
typical  path  through  the  atmosphere. 

The  determination  of  the  frequency  intervals  of  the  absorption  lines  must  be  ob¬ 
tained  in  a  laboratory  experiment  at  low  pressure  so  that  the  individual  lines  can  be 
separated.  Determinations  for  A'  =  1  to  25  by  three  investigators  are  listed  in  Table  1. 
The  results  of  the  three  measurements  can  be  considered  the  same  for  the  present  pur¬ 
pose,  and  Burkhalter’s  (6)  more  complete  results  will  be  used  in  the  computation  of  the 
attenuation. 

The  experimental  determination  of  the  constant  and  the  line  breadth  has  been 
taken  in  two  ways:  by  direct  measurement  of  each  individual  line  profile  at  low  pressui  e 
and  by  use  of  Ai-  as  a  parameter  to  obtain  the  best  theoretical  fit  to  the  measured  ab¬ 
sorption  profiles  at  atmospheric  pressure. 

The  results  obtained  for  the  line  breadth  at  low  pressure  are  consistently  about 
1.5  GHz/ atm  with  consistent  linear  pressure  dependence,  while  field  measurements  at 

Table  1 

Experimental  Determinations  of  Absorption  Line 
Frequency  Intervals  and  Line  Breadths 


Series  (a,  ) 

Absorption  Line  Frequencies, 
(GHz) 

Line  Breadth  Factors  at  Unit  Pressure, 

/'  (GHz/atm) 

Burkhalter 

Mizusuhiraa 

Gokhale 

Art  man 

Burkhalter 

Gokhale 

Anderson 

(6) 

(1) 

(10) 

(ID 

(0) 

(10) 

0) 

25- 

53.5922 

53.5994 

_ 

_ 

_ 

23- 

54.1300 

54.1294 

— 

— 

— 

— 

1.13 

21- 

54.6725 

- 

— 

— 

0.63 

— 

0.96 

19- 

55.2208 

55.2216 

— 

1.46 

— 

— 

1.23 

17- 

55.7841 

55.7846 

— 

1.38 

— 

— 

1.34 

1+ 

56.2651 

56.2656 

56.2656 

1.67 

— 

1.05 

1.49 

15- 

56.3628 

56.3640 

56.3842 

1.51 

1.44 

— 

13- 

56.9687 

56.9708 

— 

1.41 

— 

1.42 

11- 

57.6120 

57.6114 

57.6123 

1.50 

— 

1.02 

— 

9- 

58.3240 

58.3232 

58.3249 

1.52 

— 

0.63 

1.50 

3+ 

58.4462 

5S.4462 

58.4463 

1.69 

— 

0.84 

1.30 

7- 

59.1634 

59.1634 

5C.1642 

1.38 

— 

0.84 

— 

5+ 

59.610 

59.5914 

— 

1.49 

-- 

— 

— 

5- 

30.3064 

60.3080 

— 

1.51 

— 

— 

1.41 

7+ 

60.436 

60.4334 

1.46 

— 

- 

1.56 

9+ 

61.120 

61.1496 

1.47 

— 

— 

— 

11+ 

61.8002 

61.7998 

— 

— 

— 

1.50 

13+ 

62.4117 

62.4133 

62.4129 

— 

0.65 

0.75 

— 

3- 

62.4861 

62.4872 

62.4862 

— 

1.59 

1.11 

1.46 

15+ 

62.970 

62.9966 

— 

— 

— 

1.35 

17+ 

63.5683 

63.5672 

— 

— 

— 

— 

19+ 

64.1276 

64.1280 

— 

— 

— 

1.20 

21+ 

64.6789 

64.6782 

— 

— 

— 

— 

— 

23+ 

65.220 

65.2242 

— 

— 

— 

— 

0.96 

25+ 

65.770 

— 

— 

— 

— 

— 

1- 

— 

118.7505 

— 

- 

- 

- 

Pressure 

=  1  mm 

=  1 

_ 

=2  to  10 

<  1 

=1  to  6 

=  0  to  12 

(mm  Hg) 
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atmospheric  pressure  are  consistently  about  0.6  GHz/ atm  (12-15).  The  low-pressure 
individual  line  observations  are  also  summarized  in  Table  1.  This  apparent  inconsist¬ 
ency  has  been  notea  by  others  (4,16)  and  various  theories  have  been  suggested  —  notably 
tnat  the  effect  of  all  the  absorption  lines  is  not  a  simple  sum. 

Artman  and  Gordon's  data  (11)  at  intermediate  pressures  provides  a  basis  for  de¬ 
termining  and  the  actual  pressure  dependence  of  the  line  breadth  factor. 

Best-fit  values  for  ...  and  .  are  obtained  separately  for  sets  of  data  points  at  1, 
0.5,  ar.d  0.25  atm  pressure  taken  from  Artman  and  Gordon's  paper  representing  meas¬ 
urements  of  dry  air.  A  least-squares-fit  computer  program  is  applied,  assuming  that 
the  theoretical  equations  (Eqs.  (2)  and  (3))  give  at  least  the  functional  form  of  the  profile 
correctly.  The  results*  (for  7  =  300=K)are: 


"  (mm  Hg) 

(dB/km  GHz) 

■  /' (GHz/atm) 

i  760 

0.01647  t  0.00076 

0.873  t  0.195 

380 

0.01790  t  0.00153 

1.496  c  0.456 

180 

0.01742  i  0.00065 

1.518  r  0.200 

Since  .  should  be  constant,  it  can  be  fixed  at  its  mean  value  =  0.0173  and  the  best- 
fit  values  of  ' .  obtained: 


/’  (mm  Hg) 

•„  (dB/km  GHz) 

/’(GHz/atm) 

760 

0.0173 

1.0791  t  0.0585 

380 

0.0173 

1.2987  t  0.1583 

190 

0.0173 

1.4995  i  0.1659 

In  each  case  the  fit  is  of  about  the  same  precision,  and  without  further  t  -ta  it  will  be 
assumed  that  • ,  is  indeed  constant.! 

The  choice  of  .  is  quite  important,  since  the  value  largely  determines  the  absorp¬ 
tion  at  frequencies  removed  from  the  resonant  peaks,  which  is  the  region  of  greatest 
interest  for  propagation  purposes.  The  theory  of  the  line  broadening  has  been  the  subject 
of  many  studies,  but  none  has  yet  predicted  -  tactiy  the  correct  value  or  behav;or  with 
temperature  and  pressure  so  that  it  is  necessary  to  use  a  semiempirical  approach  using 
the  magnitude  of  ' .  as  a  fitting  parameter  and  a  semiempirical  model  to  give  the  pres¬ 
sure  and  temperature  dependence.  Since  the  use  of  .  as  a  curve-fitting  parameter  becomes 
quite  subjective  when  only  a  few  data  points  are  available  as  in  some  of  the  actual  path 
loss  measurements  (13,15),  the  model  for  ';</•.  n  used  in  the  present  work  is  based  on 
the  study  of  Artman  and  Gordon’s  data.  The  adopted  model  is 


♦lannurz:  (17)  claims  there  is  another  related  efiect  of  the  collision  process  in  winch  the 
phase  shifts  induced  by  the  collisions  results  in  an  effective  displacement  of  the  resonant 
frequencies  with  increasing  pressure.  He  uses  a  pressure  dependence  of  the  form 


0 


to  obtain  a  better  fit  to  measurement s  obtained  at  pressures  above  1  atm  using  a  single 
nne  at  „  -  60  GHz  to  approximate  the  entire  spectrum.  In  fitting  Artman  and  Gordon's 
data,  :  was  used  as  an  additional  parameter,  allowing  the  resonant  frequency  of  each 
resonance  to  be  shifted  by  an  amount  proportional  to  the  pressure.  The  mean  value  of  a 
obtained  for  the  six  fits  is  -0.04  ±  0.27  MHz/mm.  Without  further  evidence,  this  effect 
will  be  neglected  for  the  present  purpose  and  ;  will  be  set  equal  to  zero. 

(This  value  of  -0  is  greater  than  the  theoretical  value  of  ■„  =  0.0145,  which  i.  obtaired 
from  the  various  constants  and  an  assumed  number  fraction  o<  oxygen  (0,21). 


LINE  efiCAOTH,  Av,  <5Hl) 


TfMjp%-'t  ■»  fw.-u  i  ssn  m 
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0  00197  (.11/  Pimm  Mg).  0  <■  P  •_  $10  . 

At  -  0  690  (41/  [  /*(  mm  !!*>  ;60'1  *tJ.  / '  '  r.o  . 


I00^ - - - 

r  _ EXTRAPOLATED  LOW-PRESSURE 

t  RESULT  (IS  GHJ/ATM) 

i  _ MOOEL  USED  IN  AT  TE  VJSTIDN 


Fig.  1  -  Lin«-  breadth  vs 
pressure  modeled  toArtman 
and  Gordon's  (11)  data 


The  values  of  A:  </’>  obtained  in  fitting  Artman 
and  Gordon's  data  and  the  resulting  mouel  of  the 
pressure  dependence  of  '  are  plotted  in  Fig.  1.  It 
is  clear  that  begins  to  depart  from  a  linear  de¬ 
pendence  on  P  for  P  >  300  mm  Hg.  The  data  points 
ar  e  plotted  with  the  semiempirical  curves  that  re¬ 
sult  from  this  model  in  Fig.  2. 

The  temperature  dependence  of  .  at  constant 
pressure  has  been  measured  by  Tinkham  and 
Strandberg  (18)  (ro.  =  0.75)  and  Beringer  and 
Castle  (19)  (to,  =  0.87).  A  mean  value  of  '().  =  0.8 
wil:  be  used. 


Water  Vapor 

Available  measurements  of  water  vapor  absorp¬ 
tion  are  less  plentiful.  The  values  of  . ’ ,  Ev ,  and  7. 
for  the  resonant  water  vapor  absorption  computed  by 
Van  Vleck  (5)  are  found  to  satisfactorily  describe 
the  experimental  data  of  Becker  and  Autler: 


0. 1864 


JH  m3 
km  G1U  ft 


6  14  °K 


7.  *  170 


The  line  breadth  parameter  has  been  measured  by  Becker  and  Autler  (20)  and  is  found  to 
depend  on  the  water  vapor  density. 


=  ,A;,,o  T  n 
0 


+  0  (IP  46.  ) 


at  constant  T.  This  vapor  density  dependence  is  an  experimental  result  which  reflects 
the  greater  effectiveness  of  H20-H20  collisions  in  producing  line  broadening  compared 
to  HoO-air  collisions.  The  value  of  P  obtained  at  318°K  and  small  a  is 


For  the  temperature  dependence  of  At,  the  computed  result  of  Benedict  and  Kaplan 

(21)  will  be  used  in  the  absence  of  experimental  results:  vu  0  =  0-63. 

? 

The  residual  absorption  of  the  submillimeter  resonances  is  given  by 


p 


i/:Av  . 
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E  P  •  760  mm  Hg 
*  P  •  580  mm  Hg 


Fig.  i  -  Artman  and  Gordon's  data  (11}  with  semi- 
empirical  absorption  curves,  holding  >c  constant 
(0.0173  dB/km  GHz) 


The  value  of  ( •  ’rs>1)  is  found  to  increase  with  the  water  vapor  density.  Schulkin  (22) 
has  calculated  <  •  r  „  „  i  „  to  fit  Becker  and  Autler’s  data  (20)  using  their  result  for  the 
line  breadth: 


<> 


*  dB  m5 

(  1  83  +  0.  02,  1  -  iO'*  - 

km  g  Ohz  * 


The  water  vapor  absorption  curves  are  plotted  with  the  oxygen  curve  in  Fig.  3.  The 
curves  show  that  absorption  by  water  vapor  largely  controls  the  total  attenuation  over 
the  decade  from  10  to  100  GHz  except  the  region  of  high  attenuation  due  to  the  oxygen 
resonances. 


TOTAL  ABSORPTION  FOR  PATHS 
THROUGH  THE  ATMOSPHERE 

The  absorption  for  a  path  through  the  atmosphere  is  given  by 


.l(:.F)dIt=  Jl  . 

4) 

where  l  is  the  distance  along  the  path  measured  from  the  ground,  and  T  and  /*  vary 
with  altitude  which  in  turn  is  a  function  of  /  and  K,  the  elevation  angle  of  the  antenna. 

The  variation  of  temperature  and  pressure  with  altitude  is  not  readily  measured,  so 
that  a  model  is  needed  to  represent  typical  or  average  conditions.  The  ARDC  1939  model 
atmosphere  (23)  is  used  for  this  purpose.  A  water  vapor  distribution  proportional  to  a 
10-year  average  distribution  for  the  Washington,  D.C.,  area  (24)  is  also  assumed.  The 
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Fig.  3  -  Water  vapor  and  oxygen  absorption 
coefficients  vs  frequency  (T  -  15°C  -  bO^F. 
/*  =  7b0  rmr.  Hg,  .  0  =■  7.5  g/m  ’) 


ARDC  profiles  and  the  average  water  vapor  distribution  are  shown  in  Fig.  4.  The  alti¬ 
tude  dependence  of  the  absc’^tion  coefficient  (dB/km)  that  results  from  this  model  at¬ 
mosphere  is  r.nown  in  Fig.  F  The  curves  plotted  in  Fig.  5  illustrate  the  difference  in  the 
dependence  on  altitude  tor  frequencies  close  to  the  resonances  of  the  absorption  ve.g.,  60 
GHz)  and  those  well  removed  from  the  resonances  (e.g.,  50  GHz).  This  difference  is  due 
to  the  cancellation  of  the  pressure  dependence  of  the  absorption  for  frequencies  exactly 
equal  to  that  of  the  resonances.  For  thest  frequencies  the  absorption  coefficient  remains 
nearly  constant  to  high  altitudes,  and  for  tne  nearby  frequencies  the  decrease  in  the  ab¬ 
sorption  is  controlled  only  by  the  narrowing  of  the  line  breadth.  The  discontinuities  that 
occur  above  ~  12  km  are  due  to  the  arbitrary  cutoff  in  the  assumed  water  vapor  distri¬ 
bution. 

The  absorption  profiles  for  various  altitudes  are  three  values  of  ground-level  water 
vapor  density  (.  0),  and  are  shown  in  Figs.  6  through  8.  The  chosen  density  values  of  7.5, 
15.0,  and  22.5  g/m1  correspond  to  60-percent  relative  humidity  at  60,  80,  and  93 °F,  re¬ 
spectively.  The  curves  show  the  rapid  decrease  with  altitude  of  the  absorption  coefficient 
at  frequencies  away  from  the  resonances  and  the  nearly  altitude- independent  behavior  at 
the  oxygen  peaks.  Although  the  plotter  does  not  have  sufficient  resolution  to  reproduce 
the  full  height  of  peaks  at  the  resonances,  the  line  structure  of  the  oxygen  absorption  is 
in  evidence  at  the  highest  altitude. 

In  the  actual  calculations  applied  here,  the  atmosphere  is  civided  into  a  series  of 
110  spherical  shells  of  exponentially  increasing  thickness  to  a  height  of  30  km.  *  Except 


♦Repeating  the  calculation  wi»h  liO  layer.-  from  0  to  30  km  having  uniform  thicknesses 
produced  esseutia'ly  the  same  results. 
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Fig.  4  -  Temperature  and  pressure  profiles  of 
the  1959  AKDC  model  atmosphere,  and  average 
water  vapor  distribution 
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Fig.  5  -  Net  altitude  dependence  of  absorp¬ 
tion  for  the  1959  ARDC  model  atmosphere 


FREQUENCY  (GHz) 


Fig.  6  -  Absorption  profiles  for  various  altitudes 
for  p0  =  7.5  g/m3  (60%  RH  at  60*F) 


1000 


Fig.  7  -  Absorption  profiles  for  various  altitudes 
for  p0  =  15.0  g/m3  (60%  RH  at  80* F) 
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Fig.  8  -  Absorption  profiles  for  various  altitudes 
for  ,„  =  22.5  g/m3  (60%  RH  at  93*F» 


for  single  frequencies  corresponding  exactly  to  the  resonances,  the  contribution  of  the 
attenuation  from  altitudes  greater  than  30  km  is  negligible. 

The  slant  distance  ln  thru  the  layer  between  the  altitude  hr  and  altitude  hntl  is 
given  by 

ln  ~  ac  {Mn  sin  En 

where 


and  oP  is  the  radius  of  the  earth.  Kn  is  the  elevation  angle  of  the  ray  path  at  the  point  at 
which  it  eaters  the  nth  layer.  For  the  first  layer,  ff,  is  the  antenna  elevation  angle  and 
the  succeeding  values  of  Kn  are  given  by 


f>n  i  - 1 

a  '  ,lr.  1 1  =  a 


£,  r  cos 


n-  I 

-fC 


— il —  cos 


(F 


n  -  J 


>j 


where  \n  is  u.c  index  of  refraction  at  the  « th  altitude. 


The  integral  is  then  approximated  using  the  trapezoidal  rule 
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where  the  subscript  n  denotes  the  value  computed  for  the  «th  altitude.  Similar  calcula¬ 
tions  that  do  not  allow  for  the  bending  of  the  ray  path  show  that  the  attenuation  is  essen¬ 
tially  the  same  for  a  given  antenna  elevation  angle  and  that  the  only  effect  of  refraction 
is  to  change  the  angle  required  for  a  given  source. 

The  final  results  of  the  calculations  giving  the  total  attenuation  due  to  atmosD^rric 
oxygen  and  water  vapor  are  shown  in  Figs.  9  and  10  as  a  function  of  frequency  sor  vari¬ 
ous  antenna  elevation  angles  and  as  a  function  of  elevation  for  various  frequencies  in 
Fig.  11.  The  effect  of  the  altitude  dependence  has  sharpened  the  peaks  and  has  produced 
a  steep  increase  in  the  total  attenuation  below  10  degrees  elevation.  The  plot  in  Fig.  10 
shows  on  a  linear  scale  the  frequency  and  elevation  dependence  of  total  attenuation  (less 
than  20  dB)  and  emphasizes  the  extremely  steep  frequency  dependence  near  the  oxygen 
resonances. 


Fig.  9  -  Total  absorption  vs  frequency  for  patht 
through  the  atmosphere  at  various  antenna  ele¬ 
vation  angles  ( p0  -  7.5  g/m3 ) 


ATTENUATION  BY  RAINFALL 

The  attenuation  due  to  rainfall  is  treated  separately  since  the  distribution  of  rain 
along  the  propagation  path  is  generally  unknown  and  unpredictable,  making  the  .  asults  of 
rain  attenuation  calculations  quite  qualitative. 

The  theoretical  treatment  of  rainfall  attenuation  is  an  exercise  in  classical  electro¬ 
magnetic  theory.  The  general  solutions  to  the  wave  equation  in  expansions  of  spherical 
functions  are  known.  An  incident  plane  wave,  a  scattered  wave,  and  a  transmitted  wave 
are  expressed  in  the  spherical  functions,  and  the  coefficients  of  the  expanded  solutions 
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are  fitted  to  the  boundary  conditions  at  the  surface  of  a  spherical  raindrop.  The  attenua¬ 
tion  due  to  a  single  raindrop  is  obtained  by  comparing  the  absorbed  and  scattered  ener¬ 
gies  with  the  incident  energy. 

The  total  attenuation  due  to  the  rainfall  is  then  reiated  to  the  rainfall  rate  using  an 
empirically  based  model  of  the  drop  size  vs  number  distribution  as  a  function  of  rain 
rate  by  summing  the  attenuation  due  to  the  various  si/.e  drops,  neglecting  the  effect  of 
multiple  scattering. 

The  attenuation  at  frequencies  from  4  to  100  GHz  has  been  most  recently  treated  by 
Oguchi  (25,26)  and  has  been  reviewed  with  reference  to  the  experimental  evidence  by 
Medhurst  (27)  and  Crane  (28).  These  papers  also  contain  the  earlier  background  refer¬ 
ences,  detailed  discussions  of  the  theoretical  solutions,  and  choice  of  drop  size  vs  num¬ 
ber  vs  rate  distributions. 

In  attempting  to  estimate  the  actual  effect  of  rainfall  on  propagation  along  paths 
through  the  atmosphere,  the  gi  eatest  difficulty  is  not  in  computing  U.a  magnitude  of  the 
attenuation  once  the  characteristics  of  the  rain  are  known,  but  in  predicting  these  char¬ 
acteristics.  In  examining  the  attenuation  observed  during  a  given  storm,  the  dif  "epancy 
between  the  measured  and  expected  values  will  in  greater  part  be  due  to  departures  of 
the  rate  and  rate  vs  drop  size  distributions  from  the  average  and  to  the  vagaries  of  the 
spatial  distribution  models  in  both  altitude  and  range.  The  results  of  Oguchi's  calcula¬ 
tions  are  in  good  agreement  with  results  of  experiments  in  which  the  various  rain  param¬ 
eters  are  known,  and  these  results  may  be  used  without  adding  significantly  to  the  error 
inherent  in  the  assumption  of  models  for  rainfall  distribution.  Curves  of  attenuation  vs 
frequency  and  rate  interpolated  from  Oguchi's  tabulated  numerical  results  are  given  in 
Figs.  12  and  13. 


Fig.  12  -  Theoretical  rainfall  attenuation  vs 
frequency  for  various  rain  rates 
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Fig.  13  -  Theoretical  rainfall  attenuation  vs 
rain  rate  for  various  frequencies 


Tc  discuss  rainfall  attenuation  on  a  general  basis  without  reference  to  a  particular 
storm,  it  is  necessary  to  use  a  statistical  approach.  A  meaningful  sampling  of  typical 
rainfall  effects  may  be  obtained  by  considering  values  of  rainfall  rates  which  will  be 
exceeded  1,  0.1,  and  0.C1  percent  of  the  time  corresponding  to  light,  moderate,  and  heavy 
rain. 

Statistical  studies  of  rain  rate  as  a  function  of  the  total  time  a  specific  rate  was 
exceeded  over  the  time  period  of  a  year  are  available  for  various  locations.  The  survey 
papers  of  Weible  and  Dressel  (29)  and  Medhurst  (27)  contain  lists  of  references  for  this 
and  other  types  of  studies  that  will  be  referred  to  here. 

The  results  of  Bussey  (30)  for  Washington,  D.C.,  give  the  instantaneous  rain  rate  as 
a  function  of  the  total  number  of  hours  per  year  dur  ing  which  the  rate  is  exceeded.  The 
rate  exceeded  87.6  hours  per  year  or  1  percent  of  the  time  is  3.05  mm/hr;  15.2  mm/hr 
is  exceeded  0.1  percent  of  the  time;  and  61.0  mm/hr  0.01  percent  of  the  time. 

Altnough  it  is  known  that  light  rains  tend  to  fall  over  wide  areas  simultaneously 
while  heavy  rains  are  generally  due  to  intense,  highly  localized  storms,  studies  giving 
the  relative  rates  as  a  function  of  the  ground  range  from  the  observer  at  a  given  instant 
are  not  avaiiable.  Rather,  the  observed  quantities  are  given  in  terms  of  averages  over 
time  or  over  collecting  area.  Such  averages  necessarily  include  the  effects  of  storm 
translation  and  are  difficult  to  interpret  precisely  in  terms  of  momentary  distributions. 
Radar  studies  of  individual  storms  (31)  have  suggested  that  a  parabolic  distribution  is  a 
fair  approximation  of  that  observed: 

Kr(r)  -  1  -  tor  r  <  a 

1=0,  for  r  *  a  , 

where  /<*r(r>  is  the  rain  late  at  ground  distance  t  from  center,  relative  to  the  rate  at  the 
center,  and  a  is  the  radius  of  the  storm  (km). 


20 


R.  A.  LEFANDE 


The  dependence  of  the  rate  on  altitude  based  on  model  distributions  of  the  various 
parameters,  is  given  by  Valley  (32).  Valley's  curves  can  be  fitted  for  computational 
purposes  by  the  function 

.Va(A>  =  (  I  04$)  exp  !-0  17  4  (A  ♦  ;)-)  . 

where  ffA(A)  is  rate  at  altitude  A  (km)  relative  to  rate  at  ground. 

In  the  absence  of  meaningful  statistical  studies  of  the  instantaneous  rainfall  distri¬ 
butions,  the  following  approximate  models  will  be  assumed: 


The  cross  sections  ^  these  model  storms  are  shown  in  Fig.  14a,  b,  and  c.  To  illus¬ 
trate  the  coarseness  o'  ,  j  models,  the  cross  section  of  a  real  thunderstorm  is  repro¬ 
duced  in  Fig.  14d. 


RANGE  (km) 

a.  Model  1:  a  -  50  km 


RANGE  (km*  RANGE  (km)  RANGE  (km) 


b.  Modil  2:  a  ~  20  km  c.  Model  3:  d.  Real  thunder- 

<7  =  10  km  storm 

Fig.  14  -  Model  rainstorms  and  a  real  thunderstorm, 
showing  contours  of  rel_tivr  rain  rate 
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Since  it  is  not  possible  to  consider  possible  preferred  azimuthal  directions  without 
prior  knowledge  of  the  synoptic  conditions  at  a  given  site,  it  will  be  further  assumed,  as 
a  coarse  average  over  all  directions,  that  the  storm  is  centered  on  the  site. 

The  range  and  altitude  distributions  may  be  combined  to  give  an  effective  path  length 
through  the  storm  as  a  function  of  elevation  angle  £.  The  effective  path  length  is 


!■  </■;>  J  It, i  < ....  6' l  It, .<  :  k  K)  <H  . 

the  distance  through  an  equivalent  storm  of  constant  unit  rain  rate  throughout,  having  the 
same  total  attenuation  as  the  model  storm.  The  use  of  such  an  effective  path  length  as¬ 
sumes  that  the  attenuation  coefficient  has  a  nearly  linear  dependence  on  rain  rate.  The 
approximate  attenuation  contributed  by  the  rainfall  can  be  obtained  by  multiplying  this 
length  by  the  attenuation  coefficient  corresponding  to  the  local  rain  rate  and  the  operat¬ 
ing  frequency; 


j  K.  it ,  > 


it  i  /’  (hi 


where 

=  attenuation  due  to  rain  (dB), 

,tilB  =  attenuation  coefficient  (dB/km),  and 
V  local  raia  rate  (mm 'hr). 

Figure  15  shows  the  effective  path  length  as  a  function  of  elevation  angle  for  .  conns  of 
various  radii  based  on  the  above  models. 

The  rainfall  attenuations  computed  in  this  way  using  the  rough  models  should  indi¬ 
cate  at  least  the  order  cf  magnitude  of  the  effect  likely  to  be  encountered.  The  attenua¬ 
tion  values  for  the  three  model  storms  are  added  to  the  absorption  due  to  oxygen  and 
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Fig.  f  5  -  Effective  path  length  vs 
el'-vation  for  model  storm* 
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water  vapor  and  displayed  as  total  attenuation  in  Figs.  16  through  18.  It  is  clear  that  the 
rainfall  attenuation  is  dominant  at  the  higher  frequencies  for  any  but  the  lightest  rain 
storm. 


Fig.  16  -  Total  attenuation.  Model  Storm  1 


CONCLUSION 

The  values  of  absorption  coefficients  used  here  are  computed  in  detail  from  the  the¬ 
ory  as  tempered  by  the  available  empirical  evidence.  They  represent  as  accurately  as 
possible  the  effect  of  atmospheric  oxygen  and  water  vapor  on  microwave  radiation. 

The  values  of  total  path  attenuation,  however,  also  depend  on  artificial  models  for 
the  distribution  of  the  atmospheric  constituents,  nnd  the  results  must  be  viewed  as  rep¬ 
resentative  of  values  likely  to  be  encountered  but  not  necessarily  precise  under  all 
conditions. 

The  results  obtained  for  the  rainfall  attenuation  are  intended  only  as  a  qualitative 
indication  of  the  effect  to  be  expected.  Serious  attempts  to  predict  the  rain  losses  must 
rest  on  meaningful  statistical  studies  of  instantaneous  rainfall  distributions  at  the  site  of 
interest.  Experimental  evidence  is  limited,  and  determination  of  the  actual  attenuation 
occurring  along  atmospheric  paths  at  elevation  angles  other  than  zero  must  be  made  to 
refine  the  computed  results  and  to  determine  the  spread  of  values  under  a  range  of  oper¬ 
ating  conditions. 
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Fig.  17  -  Total  attenuation.  Model  Storm  Z 
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SYMBOLS  (in  order  of  appearance) 
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A 

J 

•“o 

^A: 

U.f.  ) 

k 

T 

Z 

/,<*» 
fr  (  J.  A,  V) 


B 

/, 


Fraction  of  molecules  in  j  th  energy  state 

Absorption  coefficient 

Frequency  of  radiation 

Number  of  molecules  per  unit  volume 

Quantum  mechanical  matrix  element  connecting  ?  th  and  /th  molecular 
energy  states 

Fraction  of  molecules  in  energy  state  < 

Resonant  frequency  of  absorption  line  corresponding  to  the  transition  from 
state  i  to  state  / 

Line  shape  function 

Rotational  quantum  number 

Total  angular  momentum  quantum  number 

The  Bohr  magneton 

Matrix  elei.  3nt  corresponding  to  transitions  A  -  A-  1  ( +  >  and  A  -  i  -  A  ( - ) 

Statistical  weight  factor 

Energy  of  ;  th  state 

Boltzmann's  constant 

Absolute  temperature 

Partition  function 

Fraction  of  states  in  vibrational  state  defined  by  principle  quantum  number 

Fraction  of  states  in  rotational  states  defined  a  set  of  quantum  numbers  •/,  / 
and  </ 

Quantum  number  for  the  projection  of  J  on  the  molecular  axis 
Base  of  natural  logai  ithnis 
A  constant 
Planck  constant 
Partial  pressure  of  o2 
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l‘  Total  pressure 

y0  Constant  adjusted  in  this  work  to  fit  the  theoretical  expression  for  the  absorp¬ 
tion  coefficient  to  experimental  values  at  pressure  i‘u ,  temperature  70 

VK'  Resonant  frequency  corresponding  to  allowed  transitions  A'  ->  A'  +  K+)  and 
A-  1  -  A'<-  ) 

x  Dummy  parameter 

^  Resonant  frequency  of  water  vapor  absorption 

<  Water  vapor  density 

*  t,  Water  vapor  density  at  /’0. 70 

Experimental  constant  fitted  to  resonant  water  vapor  absorption 
f'0  Energy  of  lower  state  of  water  vapor  resonance  transition  at  22.254  GHz 
.  Line  breadth  factor 

Mean  time  between  molecular  collisions 
Collision  cross  section 
A'  Intermolecular  force 
r  Distance  between  molecules 
n  Exponent  of  intermolecular  force  law 
( ‘.s')  0  Measured  line  breadth  at  ro .  ro 

Exponent  of  temperature  dependence  of  line  breadth  at  constant  pressure 
Jim  Atmosphere  (pressure) 

Constant  of  suggested  pressure  dependence  of  resonant  frequencies  of  absorp¬ 
tion  lines 

Value  of  -  for  o; 

'ii  0  Value  of  -  for  >L0 

•  |rs  Experimental  constant  fitted  to  residual  water  vapor  absorption 
a;  Elevation  angle  of  antenna 

l(:  .a.  )  Total  attenuation  at  frequency  --  and  elevation  angle  F. 

!  Distance  along  propagation  path 

' ,  Distance  along  propagation  path  through  the  n  th  layer 

Degrees  Fahrenheit 

<  Jit  ;>t  t  r  Kc- J vi :i 
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hn  Height  of  nth  layer 
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a  Radius  of  earth 

«’ 

yn  Total  absorption  coefficient  at  height  /<„ 

En  Local  elevation  angle  at  the  « th  layer 
Nn  Index  of  refraction  at  height  hn 
Rr(r)  Relative  rainfall  rate  at  distance  t  from  center  of  storm 
i  Relative  rainfall  rate  at  height  h 
a  Radius  of  model  rain  storm 
Pe  Effective  path  length 

■/,,  lB(v,  )  Rainfall  attenuation  coefficient  corresponding  to  frequency  rainfall  rate  H0 
R0  Local  rain  rate  (mm/hr) 

Ara in  Attenuation  due  to  rain  (dB) 
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ABSTRACT 


Microwave  radiation  is  absorbed  bjr,  atmospheric  oxygen  and  water  vapor  and  is 
absorbed  and  scattered  by  rainfall.  This  report  provides  values  of  these  losses  in 
the  frequency  decade  from  10  to  100  GHz  for  paths  through  the  atmosphere  at  eleva¬ 
tion  angles  from  0  to  90  degrees. 

The  theoretical  expressions  for  oxygen  absorption  are  assembled  in  explicit 
form.  The  resonant  frequencies  are  taken  from  published  low-pressure  measure¬ 
ments,  and  the  entire  absorption  profile  is  fitted  to  published  intermediate  and  at¬ 
mospheric  pressure  results,  with  a  constant  of  proportionality  and  the  absorption 
resonance  line  width  used  as  fitting  parameters.  The  absorption  line  breadth  is 
found  to  deviate  f.om  the  accepted  linear  pressure  dependence,  increasing  approxi¬ 
mately  as  the  square  root  of  pressure  for  pressures  above  one -half  atmospheric. 

The  results  of  similar  analyses  of  water  vapor  absorption  are  assembled  and 
combined  with  the  results  of  oxygen  absorption  analyses  to  obtain  the  complete 
absorption  coefficients  as  a  function  of  the  atmospheric  parameters. 

The  total  attenuation  due  to  absorption  in  the  dear  atmosphere  is  computed  for 
a  model  atmosphere  with  various  water  vapor  concentrations.  Results  are  provided 
in  graphical  form,  giving  attenuation  as  a  function  of  frequency  and  elevation  angle._ 
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The  prediction  of  rainfall  attenuation  is  severely  limited  by  the  absence  of 
meaningful  statistics  on  loan  rate  and  storm_structure.  Rough  estimates  for 
rain  attenuation  are  made  .1  crude  model  storms  which  occur  for  1,  0.1,  and 
0.01  percent  of  the  time.  Total  values  of  the  attenuation  for  the  model  storms 
and  the  atmospheric  absorption  are  provided. 
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1.  On  p.  2,  nine  lines  up  from  the  bottom, J:he  word 
“processing”  should  read  precessing. 

2.  On  pages  14  and  16,  the  caption  for  ^ig.  6  belongs 
to  Fig.  8,  and  vice  versa. 
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